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ABSTRACT 


A  turbvilent  film  condensation  problem  was  formulated.  The  physical 
model  and  simplifying  assumptions  are  the  same  as  those  of  Nusselt's  original 
investigation,  with  the  exception  of  including  the  turbxilent  transports.  The 
result  agrees  with  the  qualitative  results  both  of  Seban  and  Lee,  but  does  not 
agree  totally  with  Duckler. 


R.  L.  ADA^^CZAK,  Chief 
Fluid  &  Lubricant  Materials  Branch 
Nonmetallic  Materials  Division 
AF  Materials  Laboratory 
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iMTUonncTioN 

In  ooinparison  to  the  Nusselt's  original  laminar  film  condensation  theory 
(lO),  the  higher  heat  transfer  coefficient  observed  for  the  fluids  of  moderate 
and  Iiigh  Prandtl  number  and  the  lower  heat  transfer  coefficient  associated  with 
fluids  of  small  Prandtl  number  lead  us  to  consider  the  influence  of  the  turbu¬ 
lent  transports  in  condensation.  As  the  length  of  a  condenser  plate  and/or  the 
rate  of  the  heat  transfer  (condensation)  become  large,  the  thickness  of  the  con¬ 
densate  film  is  also  increased  considerably.  Consequently,  in  view  of  the  know¬ 
ledge  of  fluid  mechanics,  for  a  thicker  condensate  film  it  is  not  unreasonable 
to  conjecture  that  the  transports  of  momentum  and  Jieat  energy  are  likely  subject 
to  the  turbulent  mechanism  in  addition  to  what  is  essentially  molecular  in  nature. 

Needless  to  say,  the  actual  rational  approach  to  this  sort  of  problem 
presents  a  serious  dilemma  due  to  the  fact  that  we  have  not  yet  conquered  the 
deadlock  difficulty  of  the  turbulent  problem  in  its  own  right.  In  other  words, 
the  appearance  of  the  Reynolds  (turbulent)  shear  stress  has  been  the  source  of 
distress  in  the  field  of  turbulent  theory;  however,  the  most  practical  solution 
has  been  the  phenomenological  approach  in  which  the  analogy  to  the  random  mole¬ 
cular  motion  plays  an  essential  role  with  the  aid  of  a  judiciously  selected  eddy 
diffusivity. 

In  the  past,  few  attemps  have  been  made  in  this  direction  by  Seban  (l2), 
Rohvsenow  et  al,  (ll),  and  Buckler  (3).  Seban  and  Rohsenow  et  al.  simply  adopted 
the  well-known  piece  meal  logarithmic  velocity  profiles  of  the  bounded  flow 
for  the  velocity  distribution  within  the  condensate  film.  With  the  additional 
assumption  of  eddy  thermal  conductivity  being  equal  to  the  eddy  viscosity,  they 
carried  out  the  heat  transfer  analysis  by  a  crude  analytical  method.  In  view 


TM  ASRCN  63-33 


3 


of  the  recent  work  by  bee  (^),  it  is  now  clear  that  the  use  of  tJie  lo^farithmic 


velocity  profile  is  valid  only  in  the  case  of  a  positively  larfje  interfacial 
shear  stress.  Recently  buckler  contemplated  a  more  ambitious  program  of  actu¬ 
ally  solvim^  for  the  velocity  distribution  within  the  condensate  film  using  the 
eddy  viscosity  expressions  of  Oeissler  and  von  Karman.  However,  his  velocity 
distribution  suffers  from  the  traditional  physical  inconsistencies  of  being  dis¬ 
continuous  at  the  intersection  of  two  regions  and  nonsymmetric  (for  zero  inter¬ 
facial  shear  stress)  at  the  center  in  the  velocity  gradient.  Moreover,  the  in¬ 
terfacial  shear  stress  is  not  satisfied  correctly  for  the  non-varnishing  values, 
furthermore,  buckler  seemingly  forgot  the  salient  fact  that  the  heat  transfer 
aspect  of  the  film  condensation  problem  is  t\\ro-dimensional  in  nature,  in  the  ab¬ 
sence  of  a  similarity  transformation. 

Lee  (6)  has  recently  studied  the  solubility  of  the  equations  of  motion 
arising  from  the  Deissler  and  von  Karman  eddy  viscosities  and  concluded  that  a 
continuous  velocity  distribution,  at  least  up  to  the  first  derivative,  can  be 
obtained  v/ith  an  arbitrary  interfacial  shear  stress.  This  allows  us  to  eliminate 
the  forementionett  objections  on  the  velocity  distribution  within  the  condensate 
film  with  the  desired  interfacial  shear  stress  being  faithfully  satisfied. 

The  purnose  of  this  paper  is  to  apply  the  correct  solution  for  the  velocity 
distribution  to  a  film  condensation  problem  for  which  the  interfacial  shear  stress 
vanishes.  Obviously,  a  problem  defying  the  similarity  transformation  is  more 
involved  than  it  appears  to  be,  due  to  the  mundane  fact  that  we  must  maintain 
the  rigorosness  of  two  dimensional  formulation. 

In  retrospect,  the  general  consensus  of  the  experimenters  i.s  that  a  more 
reliaVjle  measurement  of  the  temperature  drop  across  the  condensate  film  can  be 
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nbfiervecl  when  the  thickness  of  film  becomes  large*  Certainly,  the  temperature 
dr(? )  across  the  condensate  film  is  one  of  the  most  sensitive  parameters  which 
•’•ravel, V  influence  the  reliability  of  the  condensation  heat  transfer  data,  and 
any  experiments  worthy  of  notice  should  include  such  a  measurement*  However, 
as  the  thickness  of  the  condensate  film  grows,  it  can  no  longer  be  described 
adequately  by  the  naive  Nusselt's  laminar  film  theory,  thus  presenting  the  need 
for  the  study  of  a  turbulent  film  condensation  problem* 

FORMULATIOM  OF  PROBLEtl 

Consider  a  physical  model  of  film  wise  condensation  on  a  vertical  semi¬ 
infinite  plate,  where  the  x-axis  is  parallel  and  the  y-axis  is  perpendicular 
to  the  plate,  as  shoxm  in  Figure  1*  If  the  cold  plate  is  held  at  an  uniform 
temperature  which  is  lower  than  that  of  the  saturated  pure  vapor  T^,  the 
vapor  will  form  a  condensate  film*  As  an  initial  investigation,  let  us  restrict 
ourselves  to  the  physical  model  that  was  originally  adopted  by  Nusselt*  That 
is,  assume  the  following:  (i)  consider  only  the  balance  of  acceleration  due  :,o 
gravity  and  retardation  due  to  drag  on  the  plate  in  the  equation  of  motion,  (ii> 
neglect  the  inertia  effect  and  the  interfacial  sheer  stress  at  the  outer  edge 
of  tfie  condensate  film,  (iii)  consider  the  conduction  type  mechanism  of  heat 
energy  transport  is  dominant  in  the  heat  energy  equation,  and  (iv)  neglect  the 
convective  heat  transfer  and  subcooling  in  the  condensate  film*  However,  we  do 
differ  from  the  original  Nusselt* s  formulation  by  adding  the  effect  of  a  turbu¬ 
lent  transport  in  th>’,  form  of  eddy  diffusity. 

Under  the  condition  for  which  the  above  stipulations  are  faithfully  satis¬ 
fied,  we  can  write  ine  following  equations  of  motion  and  heat  energy  within  the 
condensate  film*  0-  course,  the  usual  assumption  of  the  constancy  of  the  physi- 
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cal  properties  within  the  range  of  temperature  involved  is  evoked! 


(1) 


(2) 


iO  *x>  (u,y)) 


i  K  +  K  (u,y)) 


*Su 

sy 

zy 


g  (  S  -  y  ) 

\  ay  /o 


where  g  =  g*  (^P/  P  ),  and  V  and  yf  denote  the  eddy  viscosity  and  thermal 
conductivity  whose  explicit  functional  forms  are  yet  to  be  specified • 

In  principle,  the  above  equations  can  be  solved  for  the  velocity  and  tem¬ 
perature  distributions  in  the  condensate  film  which  are  consistent  with  follow¬ 
ing  boundary  conditions: 

u  =  0 


T  =  T 


w 


T  =  T 


at  y  =  0  and  for  x  i  0 

at  y  =  S  and  for  x  ^  0 

However,  one  immediately  notices  the  ubiquitous  appearance  of  S  in  the  system 
of  equations,  which  can  only  be  known  as  a  result  of  analysis.  In  order  to 

define  the  a  priori  undetermined  S,  we  l«t  y  ■»  I  in  equation  2  to  obtain  the 

total  heat  energy  transport  in  the  condensate  film,  that  is: 

(K  ♦K(u,y))f21).  /((  — ) 

(3)  ^y  4  '  'o 

The  growth  of  condensate  film  thickness  i  ,  which  depends  not  only  on  the 
hydrodynamics  but  on  the  phase  change  (condensation  phenomena),  can  be  described 
by  equation  3  if  one  introduced  the  following  interfacial  condition,  viz.: 


(4) 


(  K  +  H  (u,y)), 


•VdyA  c  ax  ^0 

t  p 


dy 


Substitution  of  equation  4  into  3  gives 
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(O  — Tudy  -  K  (—)  -  0 

Cp  '*0  '  dy  ^ 

liquations  1,  2,  and  5  constitute  the  system  of  equations  we  wish  to 
study.  Clearly,  they  are  reducible  to  the  original  Nusselt*s  laminar  problem, 
if  Sf  and  were  not  present.  One  can  also  show  (?)  that  the  above  set  of 
equations  could  have  been  derived  rigorously  from  the  standard  boundaiy-layer 
approximation,  if  the  foreraentioned  simnlifying  assumptions  are  properly  ap- 
pl ied . 

For  the  eddy  viscosity,  we  shall  use  the  following  expressions! 

Deissler's  ,near  the  plate  (2)! 


^  (d,y) 


2 

n  u  y 


^  1  -  ejqp  ( 


u  y  /u  )] 


0  ^  y  *  y 


and 


von  Karman’s,  away  from  the  plate  (13)! 


(7) 


i>(u,y) 


(  /  dy  )3 

(  ^2u  /  dy2  )2 


,  y*  ^  y 


vvhere  n  and  K  are  the  empirical  constants.  The  region  of  Deissler  and  von 
Karman  are  separated  at  y*  which,  in  fact,  is  the  third  disposable  constant. 
For  the  tube  flow,  Deissler  suggested  the  value  of  y*  =  /v  in 

connection  vri.th  the  asymptotic  logaritlimic  profile; however, a  slight  diff- 

* 

erent  value  of  y  should  be  used  here  so  that  the  asymptotic  logarithmic 
profile  is  preserved  for  a  large  positive  interfacial  shear  stress  (5,  6). 
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SOME  MANIPULATIONS 

Let  us  accept  the  additional  simplifyin/?;  assumption  that  the  eddy 
thermal  conductivity  is  a  constant  multiple  of  eddy  viscosity  (for  ex¬ 
ample,  unit),  then  equations  1,  2,  and  5  must  hold  for  both  regions 
referred  to  in  equations  6  and  ?.  Since  ^;(u,y)  vanishes  at  the  wall  we 
can  from  equation  1  derive  the  value  of  the  parameter  y*  =  C  //ij  (where 
C  IS  a  numerical  value),  which  demarcates  the  ranges  of  equations  6  and  ?. 
Along  the  distance  of  a  plate,  i.e.,  x-axis,  the  precise  value  of  y*  can 
not  be  assigned  ab  initio,  due  to  its  dependency  on  h  .  Nonetheless,  only 

two  cases  can  exist,  depending  on  the  value  of  $  for  a  certain  fixed  value 
of  x;  that  is,  t i  y*  and  S>y*, 

Case  I  _  ^  S  Y* 

Physically  speaking  for  this  case  the  transnort  mechanifsn  in  the  con¬ 
densate  film  is  mildly  turbulent,  i.e.,  the  so-called  sublayer  and  buffer 
regions.  The  differential  equations  for  the  velocity  and  temperature  dis¬ 
tributions  can  be  obtained  from  equations  1,  2,  and  6,  viz.s 


(B)  - 


(g/^)(i-y) 

1  ♦  (nVi)  )  ui  y  (  1  -  exp  (  -  n2uiy/»)  ))  *  ^  ^  ^ 


ax 

(9)  - 

^y 


(  /  ay 

1  +  (n2/W  )  Pj.  y  [l  -  exp(-n2uiy/  )] 


,  S  *  y* 


where  the  velocity  distribution  in  the  Leiggler* 

^ifith  the  initial  condition  of  u  =  0  at  v  = 

1 


s  region  is  denoted  by  u  . 

1 

9,  equation  8  presents  no 
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difficulty  in  the  numerical  treatment,  even  though  it  does  not  render  it¬ 
self  to  an  analytical  solution  in  a  closed  form.  For  the  temperature  distri- 
button,  equation  9  can  be  integrated  at  once  in  y  to  gives 


(10) 


T  - 


Ii  /j  1  + 


dy 


where  AT  =  T_  «  T  and  Ii  » 
s  w  X 


Jq  1  +  Pp(n2/ 1)  )uiy  [  1  -  e:q?(-n^uiy/i>  )] 

dy 


■/- 


o  1  Pj.(nV V  )u2^  [  ^  "  expC-n^Ujy/ 1>  )] 

To  obtain  the  unique  velocity  and  temperature  distributions,  we  must 
evoke  the  physical  phenomena  of  a  condensation  process.  This  can  be  done  by 
applying  the  lieibnitz  rule  to  equation  5  to  yield  an  alternate  forms 


(11) 


Cp  D  /^T 


Py  ^ 


V  ay 


The  integrand  on  the  left  hand  side  of  the  above  equation  can  now  be  obtained 
from  equation  8  v^ith  the  help  of  an  assumption  for  the  interchange bility  of 

I 

the  differential  with  respect  to  y  and  S  ,  that  iss 

(§/»>  )  {  [^  *  )niy  (l  -  exp(-n^iy/v  ))] 

(12)  -  (nVv  )(S-  y  )y  [1  ♦  ((nVv  )niy  -  1  )  exp(-n^uiy/»> )] 

[1  +  (nV'>  )nxy  (  1  -  expC-n^ujy/u  ))]  ^ 

where  =  ^u^/aS  • 

Therefore,  one  must  now  solve  equations  8,  10,  11,  and  12  simultaneously  to 
obtain  u  (x,  y),  T  (x,  y),  and  4(x). 

Case  II  —  S  >  V* 

For  this  case,  we  must  consider  two  regions  described  by  the  different 
eddy  viscosity  formulas;  namely,  equations  6  and  7»  adjoining  at  y=y*.  The 
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eountion  of  motion  valid  for  0  ^  y  ^  y*  is  identical  lo  !.h«  r.rr  vi.i  is  efjua- 
tVon  H;  however,  for  y*  6  y  4  J^the  differential  equation  must  now  be  obtained 
from  equations  1  and  7»  Let  us  denote  the  velocity  distributions 
i.n  the  Deissler  and  von  Karman  i*egions  by  and  u^,  respectively;  wa  then 
have! 


(13) 


(g/^)(S-y  ) 


^  ♦  CnVi>  )u2^y  [1  -  exp(-n2uj^y/v>  )] 


0  4  y  4  y^" 


and 

-  W/7  )(  iUj/iy  )2 

(14)  — r  ■  -r . .  —  .y‘y‘i 

^y  y(g/i>)(S  -  y  )  -  ^ugAy 


In  obtaining  equation  l4,  the  assumptions  of  ®  ^ 

were  tacitly  evoked.  Therefore,  one  must  solve  equation  13  and  14  to  obtain 
a  v.-i'ocity  distribution  within  the  condensate  film,  which  has  the  «lesiderata 
of  continuity  and  symmetry  in  velocity  gradient. 

The  temperature  distribution  can  also  be  found  by  integratiiu'  equation  2 
with  the  aid  of  equations  6  and  7»  thus? 

T  -  Tv  1  dy 

^15)  at  I2  ltPrV(u,y)A 


where  I2 


dy 

1  +  Pj.v(u,y)  /v 


,  for  which 


V  (u,  y)  takes  the  different 


forms  given  by  equations  6  and  7  depending  on  the  value  of  y. 

The  corresponding  integrand  expressions,  appearing  in  equation  11,  must  be 
derived  from  equations  13  and  14,  respectivelyJ 
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(g/v  ){  [l  +  (nVv»  )uxy  (  1  -  exp(-n^iy/w  ))] 

-  (nV»>  )(i  -  y  )  y  [  1  ♦  ((nV»^  hiy  -  D  exp(-n2uiy/D  )]  V^} 


[l  +  (nV»^  )^iy  C  1  -  expC-n^u^y/v))] 


,  0  <  y  * 


-L/of®  ft  \  fg  tZill 

■yv  l^l  V  ■  ,y  J  („  )|jy)  JUJL  iyjl 


r  ^ 


1  3/2 


,  y*  *  y  <  S 


where  and  =  ^U2/bS  •  As  before,  the  expression  for 

£(x)  can  be  derived  from  equation  11  with  the  properly  defined  functions 
connectinf^  continuously  at  y  =  y*. 

METHOD  OF  SOLUTIOM 

All  the  equations  developed  in  the  previous  section  are  too  cumbersome 
to  be  treated  analytically;  therefore,  one  must  take  a  recourse  to  the  numer¬ 
ical  method.  For  either  case,  the  temperature  function  appearing  in  equation 
11  can  be  eliminated  by  substitxiting  either  equation  10  or  15  to  yield! 


[  j^l,2  ^  ]  ^1,2  - 


(CpAT 


where  the  subscript  1  or  2  corresponds  to  cases  I  or  II,  respectively. 
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Thu'?,  in  essence,  the  mathematical  statement  of  the  problem  is  to  find 
a  solution  of  the  above  first  order  differential  equation,  i.e.,  d « /dx  = 
P(J),*nth  th.  Initial  condition  ot  0  at  x  =  0.  As  can  easily  been  seen, 
the  subtle  oonplicatlon  arises  from  the  fact  that  F(J)  must  be  found  rigor- 
ously  by  solving  the  equation  of  motion. 

From  the  numerical  standpoint,  the  following  scheme  was  adopted  for  con 
veniency:  (i)  Integrate  equation  18  to  obtain. 


(19) 


1  (Cp  AT  /a.  )  |>  X 

/  /  ^1,2  +  Ui  2(  S  )  2  dS  -  — - 

0^0  j  P 


(ii)  For  a  given  value  ofS,  the  integrand  of  the  above  can  be/ evaluated 
uniquely  from  either  equations  8  and  12,  or  equations  13,  l4,  I6,  and  17, 
depending  on  the  value  of  y*.  Having  evaluated  the  integrand  for  the  pre¬ 
assigned,  arbitrarily  small,  interval  of  S  ,  the  left  hand  side  of  equation 
19  can  then  be  computed  numerically,  (iii)  For  each  4  ,  a  corresponding  value 
of  X  can  be  found  from  the  equality,  provided  the  value  of  parameters,  (CpAT/ 
is  specified,  (iv)  However,  in  reality  we  are  seeking  the  value  of 
4  for  a  regular  interval  of  x,  instead;  This  can  easily  be  done  by  the  use 
of  an  interpolation  for  the  fixed  x  interval  with  sufficient  accuracy.  This 
sort  of  numerical  technique  seems  to  be  quite  amenable  from  the  fact  that  the 
differential  equation  is  of  separable  type# 

Let  us  now  briefly  return  to  the  discussion  of  evaluating  the  integrand 
of  equation  19  from  the  equations  of  motion.  For  case  I,  no  difficulty  arises 
on  account  of  the  regular  behavior  of  equation  8  and  12  within  the  entire 
range,  Oiyiixy*.  For  case  II,  the  velocity  distribution  must  be  described 
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by  equations  13  and  l4,  and  the  2  equations  I6  and  17. 

For  the  velocity  distribution  within  0  y  i  y*,  equation  13  can  be  solved 
numerically  with  the  initial  condition  of  u^  =  0  at  y  =  0.  However,  within 
the  range  of  y*  -  y  -  i ,  equation  14  is  expressed  into  a  corresponding  system 
of  first  order  differential  equations  as! 

aug 

-  -  P 

,  y*  ^  y  L  S 

-  (k/AT  )  p2 

I  I  ^11  m  - 

y(g/v  ){  S  -  y  )  -  P 

Usually  the  system  of  two  first  order  differential  equations  can  satisfy 
only  two  initial  (or  boundary)  conditions.  Fortunately,  due  to  the  singularity 
of  equation  21  at  y  =S being  a  nodal  point,  it  has  been  shown  that  the  follow¬ 
ing  throe  necessary  conditions  can  be  fulfilled  (5)! 

(i)  “2  ""  “1  at  y  =  y*  (continuity) 

(ti)  P  =  at  y  =  y*  (smoothness) 

(iii)  P  =  0  at  y  =  6  (symmetry) 

Thus,  the  exact  solution  is  continuous  up  to  the  first  derivative  at  y  =  y*, 
and  has  a  vanishing  gradient  at  y  =  i .  Therefore,  the  traditional  physical 
inconsistencies  on  the  velocity  distribution  are  now  removed  in  toto. 

Similarly,  the  V  must  be  evaluated  by  equations  I6  and  1?.  Again  for 

1,2 

0  i  y  y*,  the  solution  for  equation  I6  offers  no  mathematical  handicap. 


(20) 


(21) 
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On  the  other  liand  for  y*  i  y  ^  S,  equation  17  can  be  expressed  as  a  corre- 
snondingf  system  of  first  order  equations,  as: 


(22)  -  ■  Q 

ay 

»  y*  •  y  •  s 


(23) 


aQ 

ay 


){  2  [(g/y  )(  i  -  y  )  -  P]  PQ  -  p2((g/w  )  -  Qy2} 

[(g/i> )(  5  -  y  )  -  p] 


The  singular  behavior  of  equation  21  is  preserved  in  equation  23  also,  thus 
allowing  us  to  satisfy  the  following  three  boundary  conditons: 


(ii)  Q  =  av^/^y 

(iii)  Q  =  g/y 


at  y  =  y* 
at  y  =  y* 
at  y  =  S 


As  shown  in  the  Appendix,  equation  23  has  the  point  y  =Sas  a  singular  nodal 
point  so  that  the  family  of  solutions,  Q,  converges  to  g/y  as  y-»  $  , 

Having  obtained  the  solution  for  $ (x)  and  u  (x,  y),  the  temperature  dis¬ 
tribution  T  (x,  y)  can  be  computed  from  equation  10  or  15,  Therefore,  the 
solution  of  the  simultaneous  equations  describing  the  turbulent  film  condensa¬ 
tion  is  completed. 
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HRAT  Tl^A^SVV,R  R;-Sin.T3 

Even  tliouKh  tlie  velocity  and  tenperature  distributions  within  the 
condensate  {‘il'n  and  the  corresponding  film  thickness  are  the  fundamental 
functions  describing  the  condensation  phenomena,  it  has  been  the  accepted 
practice  in  the  past  to  express  it  in  terms  of  the  auxiliary  parameters, 
such  as,  k'fussel  t  and  Reynolds  numbers.  The  obvious  reason  is  that  the  actual 
measurement  of  the  steep  velocity  and  temperature  distributions  within  an 
extremely  thin  condensate  film  is  hopelessly  unyielding;  but  the  overall 
description  of  a  condensation  process  can  be  made  via  the  heat  transfer 
coefficient  and  flow  rate,  which  render  themselves  to  experimental  confir- 
mati on. 

The  usual  definition  of  condensation  heat  transfer  coefficient  is  obtained 
from  the  following  relationship: 

(?4)  H  -  k  (^T  /  ay 

Thus  from  equations  10  or  15 


This  can  be  compared  with  the  laminar  film  heat  transfer  coefficient: 

(2<^)  ho  -  k  / 

vfhar«  th«  subscript  o  refers  to  Nusselt's  laminar  problem  and 
El  ASRCN-63-33 


(27)  io  -  [l»(CpATA)  X  /Pj.  (g/i>2>jlA 


Thierefori?!,  i"?  the  ^/usselt  number  is  defined  as  Nu  =  hL/k,  the  ratio  of 
actual  Musselt  number  to  that  of  Wusselt's  laminar  case  becomes,  from  equa¬ 
tions  25  and  26: 

(28)  Nu/NUo  «  So  / 

From  the  velocity  distribution,  the  corresponding  ratio  of  flow  rates,  which 

is  defined  as  P  =  P  |  u  dy,  can  be  expressed  asJ 

o 

(.10)  T/Po  -fu  dy/(g/^)  si  /3 

Fven  though  the  expressions  like  equations  28  and  29  are  quite  instrumental 
for  the  case  of  similarity  (see  Ref.  I,  4,  7),  they  are  not  recommended  in 
our  problem  because  the  ratios  vary  along  the  x-axis  (plate  length).  Conse¬ 
quently,  the  value  of  x  must  be  treatfid  as  another  parameter,  which  would 
then  defy  the  generality  of  the  correlation  scheme. 

To  circumvent  this  inconvenience,  let  us  intro«3uce  the  following  para¬ 
meters,  Reynolds  number  and  the  average  heat  transfer  coefficient,  defined  as: 

(30)  ^®L  “  ^  ^  u  dy  I 

Ix-L  'o  I^L 

and 
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(31)  (wVg)^/\vA  -  ^/h,2 

■/itli  a  narai'olic  velocity  distribution  and  linear  temperature  distribution 
of  Nusselt's  laminar  problem,  the  following  relationship  exists; 

(33)  (  havA  -  1.U7 

DISCUSSION 

For  a  wide  range  of  values  for  P^.  and  C  aT/x  ,  u  (x,  y),  T  (x,  y), 
and  S(x)  were  computed  numerically  using  L  =  15  cm  and  V=  0,005  cm^/sec 
(a  median  value  for  water  and  liquid  metals).  In  view  of  the  previous  in¬ 
vestigation  (5),  the  appropriate  numerical  values  for  n  =  0,124,  K  =  0,4, 
and  y*  =  23  ^  /TiS  were  adopted  here  such  that  the  asymptotic  logarithmic 

velocity  profile  is  preserved  for  a  positive,  huge  interfacial  shear  stress. 

-  The  actual  turbulent  film  thickness,  4(x),  was  plotted  in  Figure  2  for 
two  widely  differing  values  of  P^,  and  the  comparison  was  made  with  the  corres¬ 
ponding  laminar  case  So(x),  For  Pj.  =  1,0,  the  turbulent  film  thickness 
deviates  from  the  laminar  case  insignificantly  because  of  S<y*;  however,  for 
Pj.  =  0,001,  the  condition  of  4>y*  prevails  for  most  of  x.  Therefore,  the 
turbulent  film  thickness  is  much  larger  than  $o(x)  due  to  the  turbulent 
mixing.  Nonetheless,  the  presence  of  turbulent  transports  contribxites  to  the 
general  increase  of  the  film  thickness  for  all  cases.  Qualitative  information 
can  be  deduced  by  simply  letting  J*»Cx^ A.  then  we  have  y*«,C“^/2  x 
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<;  has  .  numerical,  value).  J^or  =  i.o,  C  is  small  and  y*  lies  above  the 

S(x)  as  shox^  in  Figure  3  (a).  On  the  other  hand,  for  =  o.OOl  a  larger 
value  of  4(x)  (or  C)  results  in  a  small  y*  which  then  lies  below  5Cx),as  in 
Figure  3(b),  Crudely  speaking,  y*  separates  the  condensate  film  into  the 
so-called  sublayer-buffer  region  and  turbulent  core  region.  In  the  former 
re.gxon,  the  influence  of  the  Deissler's  eddy  viscosity  is  rather  unnotice- 
able;  however,  in  the  latter  region,  the  turbulent  mixing  caused  by  the  von 
Karman’s  eddy  viscosity  increases  the  film  thickness  enormously.  This  will 
in  turn  decrease  the  value  of  y*  in  the  reciprocal  relationship.  Therefore, 
it  is  more  appropriate  to  interprets  the  two  cases  of  Figure  2  as  the  differ¬ 
ence  in  the  controlling  regions,  l<y*  or  S>y^  ,rather  than  in  the  values. 

Some  typical  velocity  and  temperature  distributions  for  C  aT/x  =  o.04 
and  X  =  15  cm  are  presented  in  Figures  4  and  5.  In  Figure  4  for  Pj.  =  lo  and 
1,  the  milder  influence  of  turbulent  transports  is  manifested  by  the  nearly 
parabolic  profile;  but  as  Pr  decreases,  the  inclusion  of  the  von  Karman  re¬ 
gion  gives  rise  to  a  flatter  velocity  profile  due  to  the  violent  turbulent 
mixing.  In  contrast  to  the  previous  attempts  (3,  11,  12)*,  one  must  note 
that  the  present  analysis  yields  a  velocity  distribution  whose  gradient  not 

only  is  continuous  across  the  film  thickness  but  also  vanishes  at  the  outer 
edge  of  the  film. 

The  corresponding  temperature  distributions  are  shown  in  Figure  5,  for 


X 


:v. 


V-  ! 

f  ' 

•5 


We  shall  not  repeat  the  discussion  on  the  velocity  profile  for  a  turbulent 
film  flow  because  it  has  already  been  presented  elsewhere  (6),  in  which  the 
works  of  ref.  (3,  n,  12)  are  carefully  scrutinized. 
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which  the  general  trend  with  respect  to  Pj.  follows  what  has  already  been 
expounded  c|ua I ita lively  in  the  previous  works  (8,  12).  That  is,  as  the 
value  of  Pj.  becomes  small,  the  magnitude  of  molecular  thermal  conductivity 
is  dominant  and  moreover,  far  from  being  overshadowed  by  the  eddy  thermal 
conductivity.  Thus,  the  linearity  of  the  temperature  profile  is  preserved  in 
the  limit.  For  Pr  =  10  and  1,  the  temperature  profile  is  convex  over  the  dia¬ 
gonal  (linear  profile);  however,  as  Pj.  decreases,  the  temperature  distribution 
twists  around  the  diagonal  and  finally  coincides  closely  with  the  linear  pro¬ 
file  for  Pp  =  0.001, 

lor  the  reason  mentioned  in  the  previous  section,  we  shall  not  present 

the  condensation  result  in  terms  of  Nu/Nuo  and  r/To,  as  given  ly  equations 

28  and  29.  However,  the  correlation  of  the  average  heat  transfer  coefficient 

and  Reynolds  number  has  been  adopted  for  our  work,  and  the  Figure  6  contains 

the  result  of  computation  for  a  wide  variation  of  P...  For  P  >  a  . 

r*  r  ^  0#1,  the  aver- 

age  heat  transfer  coefficient  is  larger  than  that  of  the  laminar  case,  but  a 
decrease  is  observed  for  smaller  values  of  P^..  It  seems  apparent  that  the  aver 
age  heat  transfer  coefficient  approaches  the  asymptote  as  Pj.  becomes  vanishing¬ 
ly  small.  The  result  of  Figure  6  generally  agrees  with  the  qualitative  study 
of  Seban  (l2),  and  the  discrepancy  is  mainly  attributed  to  the  approximate 

use  of  asymptotic  logarithmic  velocity  profile  and  the  subsequent  analytical 
treatment* 

On  the  other  hand,  there  is  a  legitimate  ground  for  expecting  to  see  a 
rather  closer  agreement  with  Buckler's  no  interfacial  shear  stress  case,  be¬ 
cause  we  not  only  have  the  same  physical  model  but  also  use  identical  eddy 
viscosity  expressions.  At  most,  however,  there  can  be  a  minor  degree  of  dis- 
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ajreert^rlL  arising  from  the  difference  in  computing  the  velocity  distribution; 
namely,  hhat  ot  exact  and  approximate  methods  as  discussed  in  ref,  (6),  Figure 
7 f  containing  the  results  of  no  interfacial  shear  stress,  was  taken  from  buck¬ 
ler’s  paper,  in  which  the  curve  for  =  0,05  was  included.  In  confomity 
with  the  results  of  Seban  and  this  work,  the  average  heat  transfer  coefficient 
is  greater  than  the  laminar  case  for  Pj,  ^  1.  However,  as  Pj*  becomes  smaller, 
the  decrease  of  the  average  heat  transfer  coefficient  is  exceedingly  drastic 
so  that  there  appears  to  be  no  limit  as  Pp  approaches  zero,  A  definite  com¬ 
parison  was  made  by  carrying  out  the  computation  for  the  set  of  parameters 
used  by  buckler.  The  result  of  computation  is  shown  for  comparison  in  Figure 
7  with  n  -  0.1,  K  =  0,4,  and  y*  =  20  i>/7g?  •  Clearly,  for  Pp^l,  the  diver¬ 
gence  is  not  so  critical;  however,  for  the  smaller  values  of  P^.,  buckler’s 
curves  not  only  have  no  similarity  to  the  exact  result  but  also  decrease  with¬ 
out  limit.  Since  buckler’s  Pj.  =  0,05  curve  passes  through  the  aggregate 
of  liquid  metals  data  (9),  the  unexplainable  lowering  of  the  average  heat 
transPei*  coefficient  for  liquid  metals  condensation  was  once  considered  as 
being  caused  by  the  unaccounted  turbulent  effect.  Unfortunately,  in  light 
of  the  aV>ove  discussion,  such  a  conclusion  does  not  seem  readily  acceptable. 

The  dramatic  discrepancy  is  obviously  caused  by  the  aberrant  analysis *of 
bucklei*,  in  which  the  two  dimensional  feature  of  the  condensation  problem 
was  nob  taken  into  account  rigorously  and  the  neglect  of  the  molecular 
thermal  conductivity  in  comparison  to  the  eddy  thermal  conductivity  was 
not  justified.  Indeed,  the  latter  is  the  prime  contributor  to  the  verit¬ 
able  fact  that  the  btickler’s  heat  transfer  coefficient  is  much  smaller 
than  the  exact  result,  because  the  molecular  thermal  conductivity  is  far 
superior  to  the  eddy  thermal  conductivity  as  Pr  decreases. 
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CONCLUSIONS 

With  hope  of  olarifyinp:  the  notorious  source  of  inconsistently  small 
heat  transfer  coefficients  observed  in  liquid  metals  condensation,  the  ori¬ 
ginal  Nusselt’s  problem  was  modified  to  include  the  effect  of  turbulent  trans¬ 
ports. 

Tlie  following  conclusions  can  be  deduced? 

i)  The  thickness  of  turbulent  condensate  film  is  always  larger 
than  that  of  the  laminar  case; 

ii)  The  eddy  viscosity  tends  to  flatten  the  velocity  distribution 
in  the  turbulent  core  region,  yet  the  parabolic  profile  is  ob¬ 
served  for  the  sublayer-buffer  region; 

iii)  for  a  large  Pj.,  the  temperature  distribution  becomes  convex 
over  the  diagonal,  which  indicates  the  violent  mixing  within 
the  condensate  film;  however,  essentially  linear  distribution 
is  observed  for  a  smaller  in  which  case  the  molecular  tlier- 
mal  conductivity  is  still  dominant; 

iv)  The  heat  transfer  result  of  the  present  investigation  agree.i 
with  the  qualitative  studies  of  both  Seban  and  Lee,  but  a  con- 
suJerable  amount  of  disagreement  is  observed  with  buckler's 
work;  and 

v)  Ihe  reported  liquid  metals  condensation  data  of  Misra  and  Bonil¬ 
la  is  still  lower  than  the  present  theorgy,  in  contrast  to  buck¬ 
ler's  claim. 
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Appendix 


SINGULAR  BEHAVIOR  OF  EQUATIONS  19  AND  21 


With  the  choice  oT  the  dimensionless  coordinates  y/i  ,  y  =  P/(gJ/v>)i 

awd  f^  =  equations  19  and  21  can  be  reduced  to  the  following  dimen¬ 

sionless  forms- 


(A_l) 


11. 


-  C  f  ‘ 


yc  1  -  n  )  -  y 


(A-2) 


d»i 


-o{2  [(i-n)  -r]T<<- 1^(  i-i<  V2} 


[d-D-y] 


3/2 


where  (!  =:  K  /w  . 

The  nodal  siti^rular  beliavior  of  equation  A_i  was  investigated  previously 
by  a  graphical  method  (5).  We  shall  here  develop  an  analytical  expi'ession 
foryatil-  1,  which  can,  in  turn,  be  useful  for  the  study  of  equation  A_2  at 
the  singular  point.  The  singular  point  can  be  brought  into  the  origin  by 
letting  5=  1  in  equations  A_i  and  A_2J 


(A-3) 


d-f 

dS 


c 


JT^ 


(A-4) 


d/ 

d5 


{2  (t  -t  )  t  (  1  -  ^)/2] 
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By  1  n ing  a  naw  variably  ci)-  equation  A-3  can  be  transformed  into 

a  more  Convenient  form* 

doJ  JZj 

(A-5)  -  * - 1  (co*Oaty*0) 

Cf^ 

Per  a  .small  v^lue  of^y  ^  series  solution  of  the  above  equation  can  be 
2  L 

developed  as  cO  ~  C  +  *  *  •  ^  from  which  a  solution  of  equation  A-3  in 
an  liipliciL  form  becomes* 


(A_5) 


I  -  Tf^  +  c^'f^  ♦  ... 


As  beiove,  Ui«  nodal  sin,^ular  behavior  at  the  origin,  i.e.,'|^=  0  at 
^  =  0,  is  repeatedly  established. 

noe  near  tlie  origin,  the  solution  A_6  permits  us  to  approximate  ( J 


Sij 


-  C  andY^fi'J,  equation  A_4  can  now  be  simplified  as, 

2  0^5^  ^  )/2 


(A-?) 


5  ^ 


or. 


(A-B) 


djif  r  1  11  1 

'TT  ' 


s 


nZ 


For  a  small  valvis  ,  l/f  «  ^  ,  then  equation  A-8  can  be  simplified 

further: 


(A.9) 


dfi 


dS  2C^5 


2«  h 


(  /  -  1  ) 
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t.h<9  introduction  of  the  transformation  ^  ^  -  1,  ev-juation  A-9  reduces  to: 


U-io) 


d 


1 

2  t  ** 


The  solution  of  the  above  equation  is 5 


(A-ll)  exp  (  -  1  /  6  ) 

✓N 

whereto  represents  the  integration  constant. 

It  can  be  deduced  from  the  solution,  A_ii,that  ^  =  0  as  which  im¬ 

plies  ^  =  \  at  the  origin.  This  concludes  the  argument  that  a  solution  of 
equation  A_2  converges  to  =  1  at  >(=  1,  or  Q  =  g/i>  at  y  =  5  ,  in  unison. 

The  singular  nodal  behavior  of  equation  21  in  that  a  family  of  solutions  con¬ 
verges  to  Q  =  y  -4  can  be  established,  perhaps  by  a  an  isocline  study  as 

was'  done  for  the  Ccise  of  equation  19  (5)» 
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s 

specific  heat  at  constant  pressure 

F(  ) 

function 

r. 

g'  (  Ap/  p  ) 

g’ 

acceleration  due  to  gravity 

u 

amount  of  heat  transferred 

h 

condensation  heat  transfer  coefficient 

I 

integral  (see  equations  10  and  I5) 

K 

von  Karman's  constant 

k 

absolute  thermal  conductivity 

L 

plate  length 

n 

Deissler*s  constant 

P 

^  u  /  S  y 

Pr 

Prandtl  number  {  f K) 

Q 

»V/  ^y 

Reynolds  number  at  x  =  L 

T 

temperature 

AT 

T  ^  T 

s  w 

u 

velocity  component  in  x.-axis 

V 

'b  wfbS 

X,  y 

coordinate  system 

y* 

separation  of  Deissler  and  von  Karman  regions 
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Greeks 


r 

$ 

n 

K 

A 

S 

r 

T 

A. 

Y 

ui 

subscripts 

av 

o 

s 

w 

overscore 


flow  rate 

condensate  film  thickness 

y'h 

thermal  conductivity  (k/Cpf  ) 
laten  heat  vaporization 
absolute  viscosity 
kinematic  viscosity 
1  -  ^ 
density 

density  of  liquid  -  density  of  vapor 
shear  stress 
Q/  (r/  »>  ) 

-  1 

P/(?;i/u) 

5  -  r 


average  over  x 
Nusselt*s  laminar  case 
saturation 
wall 

turbulent  property 
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FIGURE  6 -AVERAGE  HEAT  TRANSFER  COEFFICIENT  VS  FLOW  RATE 
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FIG.  97  COMPARISON  OF  DUCKLER'S  WITH  EXACT  RESULTS 
(TAKEN  FROM  FIGURES  A  II  AND  15  OF  REFERENCE  3  ) 


